Draft version March 24, 2011 

Preprint typeset using LAT^X style cmulatcapj v. 11/10/09 



ELECTRON INJECTION BY WHISTLER WAVES IN NON-RELATIVISTIC SHOCKS 

Mario A. Riquelme 

Astronomy Department, University of California, Berkeley, CA 94720 

AND 

Anatoly Spitkovsky 

Department of Astrophysical Sciences, Princeton University, Princeton, NJ 08544 
Draft version March 24, 2011 

ABSTRACT 

Electron acceleration to non-thermal, ultra-relativistic energies (~ 10 — 100 TeV) is revealed by radio 
and X-ray observations of shocks in young supernova remnants (SNRs). The diffusive shock accelera- 
tion (DSA) mechanism is usually invoked to explain this acceleration, but the way in which electrons 
are initially energized or 'injected' into this acceleration process starting from thermal energies is an 
unresolved problem. In this paper we study the initial acceleration of electrons in non-relativistic 
shocks from first principles, using two- and three-dimensional particle-in-cell (PIC) plasma simula- 
tions. We systematically explore the space of shock parameters (the Alfvenic Mach number, Ma, the 
shock velocity, v s h, the angle between the upstream magnetic field and the shock normal, #b„, and 
the ion to electron mass ratio, nrn/m e )- We find that significant non-thermal acceleration occurs due 
to the growth of oblique whistler waves in the foot of quasi-perpendicular shocks. This acceleration 
strongly depends on using fairly large numerical mass ratios, mi/m e , which may explain why it had 
not been observed in previous PIC simulations of this problem. The obtained electron energy distri- 
butions show power law tails with spectral indices up to a ~ 3 — 4. The maximum energies of the 
accelerated particles are consistent with the electron Larmor radii being comparable to that of the 
ions, indicating potential injection into the subsequent DSA process. This injection mechanism, how- 
ever, requires the shock waves to have fairly low Alfenic Mach numbers, Ma < 20, which is consistent 
with the theoretical conditions for the growth of whistler waves in the shock foot (Ma < (mi/m e )^ 2 ). 
Thus, if the whistler mechanism is the only robust electron injection process at work in SNR shocks, 
then SNRs that display non-thermal emission must have significantly amplified upstream magnetic 
fields. Such field amplification is likely achieved by the escaping cosmic rays, so electron and proton 
acceleration in SNR shocks must be interconnected. 

Subject headings: acceleration of particles - shock waves - cosmic rays - plasmas 



1. INTRODUCTION 

Non-thermal electron acceleration is believed to be a 
universal feature of non-relativistic collisionless shocks 
both in space plasma and astrophysical environments. 
In the Earth's bow shock, for example, electrons accel- 
erated up to few tens of keV are usually observed in 
quasi-perpendicular regions, i.e, where the angle between 
the up stream magnetic field and the sho ck normal is 
> 45° (|Oka et alJI200l IGosling et al.lll989l) . Also, radio 
and X-ray observations of supernova remnants (SNRs) 
show synchrotron radiation produced by relativistic, non- 
ther mal electrons accelerated in their forw ard shocks 
(e.g. iKovama et al.lll995t IBamba eTaL|[200l [200l . 

Despite its ubiquity, the actual mechanism for shock 
acceleration of electrons is still a mystery. The 
diffusive shock acceleration (DSA) mechanism is the 

most a ccepted theo ry fo r particle acce l eration in 

shocks (lAxford et all Il977l IKrvmskvl Il977t iBell fl978l: 
iBlandford fc Ostrikerl fl978f l . This theory assumes that 
particles are scattered by MHD turbulence both in the 
upstream and downstream regions of the shock. Under 
these conditions, particles move diffusively in the shock 
vicinity, gaining energy through many crossings of the 
shock. However, in order to be able to cross the shock 
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many times, the particles need to have Larmor radii that 
are larger than the shock width, which is controlled by 
the typical ion Larmor radius. This is the still unresolved 
"injection problem" of the DSA theory, which is partic- 
ularly stringent for electrons due to their small Larmor 
radii. 

Several methods have been applied to studying elec- 
tron acceleration in non-relativistic shocks. One ap- 
proach, typically used in the context of low Alfvenic 
Mach number (Ma) shocks, uses test particle electrons 
moving in electromagnetic fields pre-determined from 
hybrid simulations (where ions are modeled as kinetic 
particles and electrons a s a massless fluid). Recently, 
IGuo fc Giacalonel (|2010D showed that electrons can be 
energized due to repeated shock crossings as they move 
along pre-existing, large-scal e magneti c fluct uations in 
perpendicular shocks. Also, iBurgess I ((2006) used the 
same method to show that shock "ripples," which are 
fluctuations on the surface of quasi-perpendicular shocks, 
can provide the scattering necessary for efficient elec- 
tron acceleration. Although in principle these mecha- 
nisms may contribute significant electron energization in 
non-relativistic shocks, they still require a population of 
electrons that are somehow injected upstream with en- 
ergies significantly above thermal. Understanding the 
origin of these particles requires self-consistent kinetic 
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calculations based on particle-in-cell (PIC) simulations. 

PIC simulations have been used by several authors 
to study the electron acceleration problem in non- 
relativistic shocks (e.g., lAmano fe Hoshino |[2007l 120091 : 
lUmeda et al. 2009). For instance, using; two -dimensional 
PIC simulations. I Amano fc Hoshino I ()2009D showed that 
efficient electron acceleration (with spectral index a = 
2-2.5) can happen in a perpendicular, Ma = 14 shock 
due to "shock surfing" of electrons on electrostatic waves 
from Buneman instability excited at the leading edge of 
the shock foot. These simulations were done for rela- 
tively low mass ratios (nii /m e = 25) and for a particular 
geometry of the magnetic field (strictly out of the plane 
of the simulation). The dependence of this mechanism 
on different shock parameters needs to be clarified, and 
in fact, in this work we find that this mechanism is not 
very efficient at realistic m ass ratios. Al s o, us ing two- 
dimensional PIC simulations lUmeda~e t al. (2009]) showed 
that, for Ma — 5 shocks, pre-acceleration in Buneman 
waves can be complemented by further energization due 
to scatteri ng at the shock r ipples, supporting the picture 
laid out by I Burgess! (|2006T > . The energy spectrum of the 
electrons in this case, however, does not correspond to a 
power law tail, but is rather described by two Maxwellian 
distributions at different temperatures. 

One of the main difficulties of PIC studies is that us- 
ing realistic physical parameters is computationally ex- 
pensive, so unrealistic approximations have to be made. 
For instance, simultaneously capturing the dynamics of 
electrons and ions requires resolving time scales as short 
as the inverse of the plasma frequency of electrons, 
uj~1, and as long as the inverse of the cyclotron fre- 
quency of ions, oj~1 (with oj p j = (^nrij e^/rrij) 1 ^ 2 and 
u> c j — ejB/m,jC, where c is the speed of light, B is 
the magnitude of the magnetic field, and rij, rrij, and 
ej are the density, mass, and electric charge of the j 
species). Given that ui p ^/(jj c ^ — (mi/me) 1 / 2 /(va/c), 
where va = B / (^imimiC 2 ) 1 / 2 corresponds to the Alfven 
velocity of the plasma, the computing time in PIC simu- 
lations is usually reduced by using artificial mass ratios 
m.i/m e >C 1836 and/or relatively large values of va/c. 
Also, these studies are usually made using one- and two- 
dimensional simulations. Although much of the relevant 
physics can still be revealed using these approximations, 
a complete understanding of the role played by the cho- 
sen parameters is crucial before any extrapolation to re- 
alistic setups is made. Thus, in this work we present 
a systematic exploration of the space of shock parame- 
ters, paying special attention to the role of artificial mass 
ratios mi / m e on the possible electron energization. 

The effect of using a reduced ion to electron mass 
ratio mi/m e has already been highlighted by previ- 
ous PIC studies of qua si -perpendicular shock struc- 
ture dScholer et al.1 [2001 iScholer fc Matsukivo I 120041 : 
Hcll inger et al.ll2007fl . These works have shown that ex- 
cessively small mi/m e may suppress the appearance of 
oblique whistler waves in the foot of low Ma shocks, 
which would lead to periodic shock reformation. 

In terms o f ele ct ron energization, the work s of 
IScholer et all (pool : IScholer fc Matsukivo I (|200l also 
show that whistler waves can lead to significant electron 
heating in the foot of the shocks, but the investigation 
of their possible role in non-thermal electron injection 



remains to be realized^. 

In this work we study the injection of non-thermal 
electrons in quasi-perpendicular shocks, using two- and 
three-dimensional PIC simulations. We systematically 
test different regimes for the shock velocity, v s h, the 
plasma magnetization (quantified in Ma), the ion to elec- 
tron mass ratio rrii/m e , and the angle between the shock 
normal and the ambient magnetic field, 0Bn- We show 
that the presence of oblique whistler waves in the shock 
foot does lead to significant electron acceleration with 
spectral index a ~ 3 — 4. This acceleration appears to 
prefer small values of the M J 4/(m i /m e ) 1 / 2 ratio, and re- 
quires 0Bn 7^ 90°. These two conditions likely explain 
why this effect was not seen by previous PIC studies of 
this problem, where strictly perpendicular shocks with 
fairly low mass ratios were used. For the realistic value 
rrii/m e — 1836, this acceleration would require low Mach 
numbers, Ma < 20, in order to explain the electron in- 
jection fraction inferred from broadband observations of 
SNRs. Thus, if this mechanism happens to be the only 
possible solution for electron injection into the DSA pro- 
cess in SNR shocks, the Ma < 20 condition would imply 
strong amplification of the upstream magnetic field of 
these shocks. 

This paper is organized as follows. In Sj2]we explain the 
numerical setup of our simulations. In we describe the 
mechanism that brings electrons to non-thermal energies, 
by analyzing the evolution of the accelerated particles in 
one of our simulations. In []4l we explore different shock 
parameters, and determine the regimes where electrons 
acceleration occurs. Our final discussion and conclusions 
are presented in [J5l 

2. SIMULATION SETUP 

We use the el e ctromagnetic PIC code TPJSTAN-MP 
(|Bunemanl[l993l : iSpitkovskvl 120051 ) in two and three di- 
mension^. A shock wave is produced by reflecting a cold 
high-speed electron-ion plasma off a conducting wall. 
The shock forms due to the interaction between the in- 
coming and reflected beams and propagates away from 
the wall along the ir— direction. The incoming plasma 

carries a uniform magnetic field, Bq, forming an angle 
&Bn with the shock normal (which coincides with the 
x— axis). The shocked plasma stays at rest with respect 
to the box, so the simulation is done in the downstream 
frame. In two dimensions the simulation box consists 
of a rectangle in the xy plane with periodic boundary 
conditions in the y— direction. 

The incoming plasma is injected through a "moving in- 
jector" , which recedes from the wall in the x— direction 
at the speed of light. The simulation box is expanded 
in the +x— direction as the injector approaches the right 
boundary. This way the memory and computing time are 

1 We note that ILevi nson (1992, 1994) discussed an analytical 
theory of electron injection due to whistler waves excited by the 
returning electrons. These works predicted that electrons should 
be accelerated in shocks with large Mach numbers, Ma > 43. How- 
cvcr, both our simulations and the recent work by Kato & Takabc 
(2010) show neither whistler excitation nor significant electron ac- 
celeration in such shocks, raising questions about the consistency 
of the assumptions that went into the theory. 

2 In both two- and three-dimensional simulations the three 
components of particle velocities and electro-magnetic fields are 
tracked. 
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saved, while following the evolution of all the upstream 
regions that are causally connected with the shock. Fur- 
ther numerical optimization can be achieved by allow- 
ing the moving injector to periodically jump backward 
(i.e., in the — x— direction), rese tting the fields to its right 
(see lSironi fc Spitkovskv Il2009f ). Since the shock travels 
much slower than the speed of light, without a jump- 
ing injector (i.e., with the injector only moving at c) the 
upstream region would comprise most of the simulation 
domain. However, since we expect the electron accel- 
eration to occur on scales close to the shock foot (i.e., 
at a distance comparable to the typical Larmor radius 
of ions, Rh,i), keeping an upstream size of a few shock 
foot lengths should be enough to capture the relevant ac- 
celeration physics. Thus, in our simulations the moving 
injector jumps backward every 2000 time steps, with the 
first jump happening after 4000-12000 time steps (de- 
pending on the length of the shock foot). The injector 
jumps backward a distance such that its average veloc- 
ity is close to the shock speed. This method also helps 
to reduce particle heating by the numerical cold beam 
instability, which can happen after a long distance prop- 
agation of the cold plasma over the numerical grid. 

We ran a series of two- and three-dimensional simula- 
tions to explore different shock regimes. The run parame- 
ters are summarized in Tables [1] [21 and|3J The numerical 
value of the speed of light is set to 0.45 cells/time step 
in all the runs. In the two-dimensional runs presented in 
the main part of the paper, Bq is always in the plane of 
the simulation (xy plane). This condition is changed in 
Appendix [A] where the case of Bo quasi-perpendicular 
to the simulation plane is considered. 

3. THE ACCELERATION MECHANISM 

In this section we disentangle the process that gives rise 
to the electron energization, focusing on one of our shock 
simulations where significant acceleration is observed. 
We use the two-dimensional run 2D-3, with v s h/c = 0.14 
(as seen from the upstream medium), va/c = 1/50, 
mi/m e — 400, Ma — 7, and Obu = 75° (the other pa- 
rameters are specified in Table [TJ . First we will describe 
the basic features of the shock, and then we will focus on 
the process of electron energization. 

3.1. Shock Structure 

Quasi-perpendicular shocks are characterized by the 
presence of the so-called shock foot. The foot is defined 
by the existence of a beam of ions reflected by the shock, 
whose bulk velocity in the upstream frame is close to 
v s h and has comparable x and z components. The 
foot region can be seen in panels a) and b) of Figure [TJ 
(between x = 410 and 470 c/u> Pye ), which depicts the 
phase space distribution functions of the ions. This 
region covers a distance comparable to the ion Larmor 
radius, Rh,i (calculated with the upstream magnetic 
field Bo), and is located right in front of the shock 
density jump or "ramp" , as shown in panel h) of the 
same Figure. An important feature of the shock foot is 
the presence of oblique electromagnetic waves that grow 
on scales of ~ 10 c/uj Pte , forming an angle of ~ 45° with 
the shock normal. These waves can be seen in panels 
e) — g) of Figure [JJ (between x = 410 and 470 c/uj p , e ), 
which show the magnetic field along the x— direction, 
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Fig. 1. — The shock structure for run 2D-3 (Table[TJ at tuj p , e = 
10000 (tU3 c i = 10, where 0} o i is defined in the upstream), whose 
basic parameters are u a j, = 0.14c, Ma = 7, rrii/m e = 400, and 
@Bn = 75°. The phase space distribution functions of the ions, 
on the x — p x and x — p z planes are shown in panels a) and b), 
respectively. The corresponding phase space distribution functions 
for the electrons, f e , are depicted in panels c) and d), respectively. 
The distribution functions are normalized by their maximum value. 
Panels e), /), and g) show the magnetic field components along the 
x— and z— axes, and the ion density, respectively. Panel h) shows 
the one-dimensional ion density profile averaged on the y axis. 
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B x , the magnetic field along the z— direction, B Zl and 
the ion density, Uj, respectively. The oblique modes 
have a right-handed circular polarization, and their 
phase velocity is comparable to the speed of the shock. 
All these features allow them to be identified as electron 
whistler waves propagating obliquely with respect to the 
background magnetic field. 

The whistler waves have been studied in the con- 
text of shock structure evolution, using one-, two- 
, and three-dimensional PIC simulations, as well as 
hybrid codes. The two-dimensional studies show 
that whistler waves can play an important role 
by suppressing the self-reformat ion of perpendicular 
and quasi-perpendic u lar shocks (iHellinger et al.l 120071 : 
ILembege et all [200l lYuan et alj I2009Q . This effect, 
however, would be less important in the fully three- 
dimensional geometry, as shown by a recent PIC study 
(|Shinohara et al. 11201 If ). 

The exact generation mechanism of these whistlers 
is still subject to debate. One candidate mechanism is 
the so-called modified two-stream instability (MTSI), 
driven by the relative cross-field velocity between 
the electrons an d the ions in the foot of the shock s 
(|Wu et all [19831 IMatsukivo fc Scholer I [20031 [2006h . 
Since in quasi-perpendicular shocks a significant frac- 
tion of the ions are reflected into the upstream, the 
MTSI can be driven by the relative motion between 
electrons and either the incoming or reflected ions. 
These two possibilities are usually referred to as MTSI1 
and MTSI2, respectively. The analytic dispersion 
relation calculations show that the MTSI1 and MTSI2 
will grow if cos(#) > 4(1 — r)M J 4/(m i /m e ) 1 / 2 and 
cos(6>) > 4rAf J 4/(mi/rn e ) 1 / 2 , respectively, where 9 is the 
angle between the magnetic field and the wave vector of 
the waves, and r is the fr action of reflected ions, wit h 
typical values of r ~ 0.2 ([Matsukivo fc Scholer I l200l . 
Independently of which of these two varieties of the 
MTSI is more dominant, these conditions show that the 
smaller is the ratio Ma/ (mi/rrie) 1 / 2 the larger is the 
range of 9 where the excitation of whistler waves in the 
foot of shocks would be possible. 

Another possibility for whistler generation is that these 
waves can be an intrinsic compon ent of oblique quasi- 
perpe ndicular shocks. Indeed, iKrasno sclskik h et al.1 
(2002) proposed an analytical model for the structure 
of these shocks where the shock ramp is treated as 
a nonlinear whistler wave. This m odel shows that if 
Ma < M w = cos(9 B n)\(\ / m i /m e ) 1 / 2 /2, where 9 Bn is 
the angle between the shock normal and the upstream 
magnetic field, Bq, the shock precursor would contain 
a stationary whistler wave train. If M w is exceeded, 
these whistler waves would become non-linear and 
would be rather confined to the shock ramp. But if 
M A > M nw = | cos(9 Bn )\( v / m l /m e ) 1 / 2 /2 1 / 2 a station- 
ary solution for the shock structure would no longer 
be possible, which would set the condition for shock 
reformation. In particular for Run 2D-3 (9 ps 30° and 
9 B n = 75°) the condition for the whistler generation 
due to the MTSI (in any of its two varieties) appears 
to be less stringent than for the model proposed by 
IKrasn osclskik h et al.1 (|2002). This is in general the case 
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Fig. 2. — The energy spectra at different positions for electrons 
and ions (solid and dashed lines, respectively) are shown for the 
shock transition region of run 2D-3 at tui C} i = 10. The spectra are 
measured at the x/c/u) Vt< , = 330, 370, 430, and 460 positions, shown 
in black, red, green, and blue lines, respectively. These positions 
are also marked by vertical lines in panel h) of Figure [T] whose 
colors match the ones used for the corresponding spectra. 

for the simulations presented in this paper. But, even 
for runs where cos(0s„) > cos(#), the growth of whistler 
waves in the foot of quasi-perpendicular shocks appears 
to be favored for small values of the Ma/ '(m-i/me) 1 / 2 
ratio. 



3.2. Electron Spectrum 

The phase space distribution functions for electrons is 
depicted in panels c) and d) of Figure [1] Although the 
electrons are mainly heated at the shock, significant elec- 
tron energization also occurs in the foot region. This is 
seen from the electron spectra shown in Figure [21 which 
are measured at several positions in the shock region. 
The positions are marked by the vertical lines in Figure 
[TJi, with colors matching the ones of the corresponding 
spectra in Figure [2j We see that the two energy dis- 
tributions measured in the downstream (red and black 
lines) show a high-energy power law tail with spectral 
index a ~ 3.6 (for comparison a Maxwellian distribution 
corresponding to the downstream electron temperature 
is shown in black dotted line). The two spectra mea- 
sured in the foot and the beginning of the ramp (blue and 
green lines, respectively) do not show a power law tail, 
but a bump in the high energy part of the distribution. 
This high energy bump is also observed further upstream, 
showing that some accelerated electrons also outrun the 
shock as they move along the upstream magnetic field. 
The maximum Lorentz factor of the electrons (r e ~ 50) 
is consistent with the electron Larmor radii being compa- 
rable to that of the ions, i.e., T e v s h(mi/m e )/c. Also, 
two downstream ion spectra, measured at the locations 
of the vertical black and red lines in panel h) of Fig- 
ure [U are depicted by dashed lines. These spectra show 
an incomplete ion thermalization, and essentially no ion 
acceleration, with only a small, exponentially decreas- 
ing bump at energies a few times above thermal. Notice 
that the downstream electron temperature is close to the 
mean kinetic energy of the ions (for mi/m e = 100), im- 
plying energy equilibration between both species. Fi- 
nally, the spectra depicted in Figure [2] are fully evolved, 
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in the sense that they do not show significant variations 
in time. This will be the case in all the simulations shown 
in this paper after a time of ~ 6w~^ . 

3.3. The Physics of Acceleration 

In this section we illustrate the physics of the electron 
acceleration process by focusing on the evolution of a 
typical non-thermal electron. In the three panels of Fig- 
ure El we plot the energy of the electron, represented by 
the solid, black line. In order to identify the source of the 
energy gain, we also plot the accumulated energy gain, 
ej , due to the work done by the electric field along the j 
axis (red, green, and blue lines represent j = x, y, and 
z, respectively). 

This energy gain is calculated as 6j = J v e jEje/m e dt, 
where v e .j and Ej correspond to the electron velocity and 
electric field along the j axis, respectively. Since these 
cumulative energies are shown on the logarithmic scale, 
we plot their absolute values, and use dotted lines when 
they correspond to negative quantities. In each panel 
we also highlight a small time interval, whose starting 
and final points are marked by yellow and red circles, 
respectively. The particle trajectories corresponding to 
each of these intervals are shown on top of Figures [S] and 
[6] by a solid, black line. The black circle on top of the 
highlighted intervals marks the time of each field snap- 
shot. The initial increase in the electron energy happens 
in the first highlighted interval, shown in Figure Eh, with 
the particle location depicted in Figures Eh and Eh- We 
can see that this initial energization happens in the shock 
foot as the particle moves through the region of ampli- 
fied whistler waves. The energy gain is initially due to 
increase in e y , with e x contributing most of the energy 
gain after tuj p ^ e ps 1960. Although before the initial en- 
ergy growth (the one dominated by e y ) e x has already 
increased significantly, its contribution to the total en- 
ergy of the electron is almost exactly compensated by a 
decrease in e z (see rising dotted blue and solid red lines 
between tw p<e = 1737 and 1920 in Figure Eh) . Indeed, 
the particle motion until this point was dominated by 
the E x B drift, which does not allow any net work to be 
done by the electric field. However, as soon as electric 
field fluctuations on scales comparable to the electron 
Larmor radius appear, this almost perfect cancellation 
stops, and a net energy gain becomes possible. Since the 
y— axis is quasi-parallel to the initial field .Bo, the initial 
energization due to E y already suggests an important fea- 
ture of the whistlers waves: their electric field has a non- 
negligible component along the magnetic field direction. 
This parallel electric field component can be explained by 
the obliquity of the waves with respect to the initial mag- 
netic field, which makes the projection of whistler electric 
field on Bo have a magnitude ~ sin(0)|.E|. The energy 
gain due to the parallel electric field can be verified from 
Figure 2J where we have plotted in red and green the en- 
ergy gain due to the electric field perpendicular and par- 
allel to B (ej_ and en, respectively). Indeed, we can see 
that en is what dominates the initial energy gain of the 
electron, confirming that the presence of the electric field 
parallel to the magnetic field is the essential component 
in the initial electron energization. This feature makes 
this process fundamentally different from fast Fermi ac- 
celeration where the shock acts as a magnetic mirror, en- 



ergizing the electrons due to the electric field component 
perpendicular to the magnetic field. Finally, we note that 
this initial energization is described in the frame where 
the downstream plasma is at rest (usually called normal 
incidence frame). The picture, however, should be essen- 
tially the same if described from any other frame like, for 
instance, the de Hoffman- Teller frame. In particular, the 
initial energization will continue to be dominated by the 
electric field component parallel to the magnetic field. 
This is due to the relativistic invariance of E ■ B, which 
would make E ■ B/\B\ a nearly invariant quantity under 
any non-relativistic frame transformation. 

After the initial energization driven by the whistlers, 
the acceleration due to E z (the convective electric 
field) starts to play a more relevant role, as can be 
seen from Figure Eh We see that in the highlighted 
interval (marked by the two vertical, dashed lines) the 
energy increase is driven by a jump both in e z and 
£ y . In Figures Eh and Et> we can see that during these 
energy increases the electron moves mainly along the 
+y— direction, and stays close to the shock ramp. The 
relevance of the energy gain due to E z is more obvious at 
later times, when the electron Larmor radius has become 
comparable to a sizable fraction of the foot length. As 
shown in Figure Eh most of the energy is gained due to 
an increase in e z . The particle trajectory corresponding 
to the marked interval (between the two vertical, dashed 
lines) is plotted in Figures Eh and Eh- We can see 
that, as it moves (mainly along the +y— direction), the 
particle jumps a couple of times from the upstream to 
the downstream and vice versa, gaining energy in a way 
similar to shock drift acceleration. 

Apart from dominating the initial electron energiza- 
tion, the energy gain along the magnetic field, e \ , also in- 
creases the particle's ability to move along the magnetic 
field lines. This is important because if an electron veloc- 
ity along Bq becomes larger than v s h/ cos(#b„), it will be 
possible for the particle to move along the +x— direction 
faster than the shock. When this happens, the electron 
will be able to escape upstream, increasing its chances to 
be further scattered by the whistlers waves. This quali- 
tative picture already shows that having 6bu 7^ 90° is an 
important component in the acceleration process. This 
angle dependence of the acceleration is confirmed in the 
following section where the importance of the shock pa- 
rameters: 0Bn, mi/m e , v s h, and Ma for electron accel- 
eration is studied. Finally, it is important to point out 
that the evolution of thermal electrons is fundamentally 
different to the one of the non-thermal particles, which is 
what we described in this section. Thermal electrons gain 
some energy due to a random combination of the paral- 
lel and perpendicular electric field mainly in the ramp 
and shock transition region, and then are rapidly incor- 
porated into the downstream medium. 

4. ELECTRON ACCELERATION REGIME 

In this section, we seek to determine the physical 
conditions under which the electron acceleration due to 
whistler waves happens. To do this, we analyze the ac- 
celeration behavior by varying different shock parameters 
one by one. 
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Fig. 3. — The time evolution of the energy of an accelerated electron, T e — 1 (where T e is the electron Lorentz factor), is depicted by 
the solid, black line in panels a), ft), and c), which show different time intervals. Besides the electron energy, the accumulated energy 
gain, €j, due to the work done by the electric field along the j axis is also plotted (red, green, and blue lines represent j = x, y, and z, 
respectively). Since these cumulative energies are shown on the logarithmic scale, only their magnitudes are plotted. Thus, in order to 
keep the sign information, the negative values of tj < are plotted using dotted lines. The vertical, dashed lines in panels a), ft), and c) 
mark the intervals corresponding to the electron trajectories tracked in panels a), ft), and c) of Figures [5] and [5] 
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Fig. 4. — The time evolution of the energy V e — 1 (where T e is 
the e lectron Lorentz factor) of the accelerated electron analyzed in 
£|3.3I for the same interval shown in panel a) of Figure [3] is plotted 
in black solid line. Besides the total energy, the cumulative energy 
gain due to the work done by the electric field perpendicular and 
parallel to B (ej_ and en, respectively) is also plotted in red and 
green lines, respectively. 

4.1. Angle Dependence 

We illustrate the importance of the angle Obu for elec- 
tron acceleration using Figure [7] The downstream spec- 
trum of the electrons at tu} p , e = 10000 (ta> c ,i = 10) is 
plotted for simulations that only differ in their angle Osn- 
The rest of the parameters are the same as for simula- 
tion 2D-3 in Table [TJ We can see that the hardest spec- 
trum is reached for 0B n = 70° (red line), with a ~ 3.6 
while the softest spectrum is obtained for the 9bu = 90° 
case (purple line), which corresponds to a Maxwellian 
distribution with a small increase in the number of par- 
ticles at energies a few times above thermal. This re- 
sult confirms the picture suggested in §3.3[ in which the 
electron energization strongly depends on the ability of 
these particles to propagate in front of the shock as they 
slide along the magnetic field lines. This condition im- 
plies that the most accelerated particles must move pref- 



erentially along the +y— axis. This is in fact what we 
noticed when tracking the trajectory of the accelerated 
electron analyzed in ^3.31 which represents the typical 
behavior of the most energetic particles in our runs. For 
Sbu < 70° the softer spectra can be explained by the 
less efficient confinement of the particles to the shock 
vicinity. A quasi-perpendicular configuration keeps the 
electrons close to the shock, not allowing them to easily 
escape upstream. From that perspective, it is reasonable 
to think that the most efficient acceleration will happen 
when cos(#b„) is just large enough to allow the electrons 
to move at v s h in the +x— direction, with decreasing 
efficiency for larger values of cos(#b„). 

4.2. Mass Ratio Dependence 

We now explore the ion to electron mass ratio mi/m e 
dependence of the electron acceleration. To do this, we 
use three runs in two dimensions with the same param- 
eters Ma — 7, Osn — 75°, and v s h = 0.14c, but with 
mi/m e = 25, 100 and 400 (run 2D-7, 2D-8, and 2D- 
3 of Table [TJ respectively) . Their downstream spectra 
are shown by the solid lines in Figures [8jx-[8j:. We can 
clearly see the role played by mi/m e in the electron ac- 
celeration. While for the mass ratios m,i/m e = 25 and 
100 cases, the spectral index a ps 5.4 and sa 4.7, re- 
spectively, the rrii/m e — 400 case shows the decrease 
of the spectral index to a rj 3.6 0. In order to show 
the correlation between this acceleration and the growth 
of the whistler modes, in Figure [9] we have plotted the 
electric field along x, E X) in the shock transition re- 
gion of these three simulations (panels a), b), and c) for 
mi/m e =25, 100, and 400, respectively). These plots 
show how the shock foot becomes dominated by large 
amplitude whistler waves as the mass ratio increases from 
mi/m e = 25 to 400. Also this Figure confirms the length 

3 Since the typical time scale for quasi-perpendicular shock evo- 
lution is given by the ion cyclotron period, , the comparison 
between simulations of different mass ratio is performed at a fixed 
time tu c % = 10. We use this criterion because, for all the mass 
ratios in our study, tuj c ^ = 10 corresponds to the typical time at 
which the downstream spectrum becomes homogeneous, with no 
substantial variations as a function of position. 



Electron Injection in Non-relativistic Shocks 



7 



1 10 



E x /E at t« =1950 




24.9 I 



0.0 



-24.9 I 



110 




40 80 120 160 

E x /E ot tcj pie =3370 



1 10 



50 90 130 

x /(c/« p ,e) 



E x /E at t« =4840 



1 70 




I 2.0 



0.0 



1 2.0 I 



13.0 



0.0 



1 3.0 I 



120 



160 



200 

x /( C /"p,e) 



240 



280 



Fig. 5. — The electric field along x, E x , for the two-dimensional 
run 2D-3 (M A = 7, v ah = 0.14c, 8 Bn = 75°, and mi/m e = 400). 
Panels a), 6), and c) correspond to times tui v>£ = 1950, 3370, and 
4840. The field is normalized in terms of Eq = BoVi n /c, where 
Vi n is the speed at which the upstream plasma is injected, as seen 
from the downstream medium (for run 2D-3, Vi n = 0.1c). On top 
of each panel, the trajectory of an accelerated electron is tracked 
by the black line, with the yellow and red circles marking the initial 
and final time of each trajectory. The black circle, marks the time 
of each field snapshot. 



scale of the whistlers waves (~ 10c/uj Pte ) determined in 
£j3j which is consistent with previous dispersion relation 
calculations ()Wu et al.lll983[ ). In this study of mass ra- 
tio dependence we also want to include the possibility of 
electron energization by "shock surfing" of electrons, due 
to electrostatic waves produce d by the Buneman instabi l- 
ity in the foot of the shocks (jAmano fc Hoshino 1 12009). 
These waves appear in the leading edge of the shock foot 
on scales of ~ c/ui Pte , and are due to the relative velocity 
between th e upstream elect rons and the ions reflected by 
the shock (|Bunemanl 119581) . As electrons get scattered 
by these waves, they gain energy due to the work per- 
formed by the fluctuating electric field, combined with 
the convective field of the upstream plasma. Since the 
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Fig. 6. — Same as in Figure[5]but for the ion density, n^. The 
density is normalized in terms of the upstream ion density, n^o- 
The arrows represent the orientation of the magnetic field on the 
xy plane. 

wave vector of the fastest growing Buneman mode, kb un , 
is parallel to the ion beam, the growth of the waves and 
their effect on the electron acceleration are better re- 
solved when the ion motion is parallel to the plane of 
the simulation (xy plane). This is achieved if the mag- 
netic field, Bq, is quasi-perpendicular to the xy plane, 
implying that the Buneman acceleration will be sup- 
pressed in our two-dimensio nal runs with the magneti c 
field in the simulation plane (jAmano fe Hoshino 1120091 ). 
We used this behavior to distinguish the contribution of 
the "shock surfing" mechanism to electron acceleration 
by comparing the results of the two-dimensional simula- 
tions with analogous three-dimensional runs, where the 
kbun is resolved by adding a third dimension of a few 
c/ijJ P ,e- The corresponding downstream spectra of the 
three-dimensional runs with mi/m e = 25, 100, and 400 
(runs 3D-1, 3D-2, and 3D-3 of Tabled respectively) are 
depicted using dashed lines in Figures HJz, [8f>, and [8^. 
We can see significant differences between the two- and 
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Fig. 7. — The downstream spectrum of the electrons at tu) p , e = 
10000 (tw Cl i = 10) is plotted for simulations like run 2D-3 (M A = 7, 
v sh = 0.14c, 9g n = 75°, and mi/m e = 400), but using different 
angle 0g n . The purple, orange, red, blue, and black lines represent 
the cases 9 Bn = 90 (run 2D-1), 80 (run 2D-2), 70 (run 2D-4), 
60 (run 2D-5), and 45 (run 2D-6). For comparison a Maxwellian 
distribution is shown in black dashed line. 

three-dimensional spectra for rrii/m e — 25, implying an 
important contribution to acceleration due to the shock 
surfing mechanism. Indeed, when analyzing the trajec- 
tories of the energetic particles of the three-dimensional 
run with mi/m e = 25, we find that most of their energy is 
gained as the electrons get scattered by Buneman waves 
in the leading edge of the foot. The difference between 
two- and three- dimensional runs, however, is substan- 
tially reduced for the case rrii/m e = 100, and almost 
disappears for m,i/m e = 400, showing that the relative 
importance of Buneman acceleration, compared to the 
energization due to whistler waves, decreases for more 
realistic mass ratios. Further details on this mass ratio 
dependence of the shock surfing mechanism are given in 
Appendix [3] where the acceleration of electrons is stud- 
ied for two-dimensional runs with B$ quasi-perpendicular 
to the plane of the simulation. The results presented 
in Appendix [3] also demonstrate that the acceleration 
due to whistler waves disappears if the direction along 
Bo is not resolved by the simulation, which makes the 
two-dimensional configuration used in the main part of 
this paper (i.e., with Bq in plane), the most suitable for 
the study the electron injection due to whistlers. Also, 
the convergence between the two- and three-dimensional 
runs depicted in Figure [5J: shows that whistler accelera- 
tion is well modeled by our two-dimensional simulations, 
with no appreciable changes due to the additional third 
dimension of a few c/ui Pie . This result, however, does 
not rule out significant three dimensional effects when a 
scale comparable to a few whistler scales (i.e., a few 10 
c/u) Pte ) is resolved along the third dimension. We leave 
this possibility as a subject of future investigation. 

4.3. Shock Velocity Dependence 

It is interesting to see how this electron acceleration 
may change if a different shock velocity is used. We test 
the shock velocity dependence by running simulations 
with the same parameters as for the two-dimensional sim- 
ulations shown in Figure[8]but using v s h = 0.042c, which 
corresponds to Uj„ = 0.03c (simulations 2D-18, 2D-19, 
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Fig. 8. — The downstream electron spectra for simulations with 
Ma = 7, 0Bn = 75°, and = 0.14c, but different ion to elec- 
tron mass ratios, m,i/m e = 25, 100, and 400, are depicted in 
plots a), b), and c), respectively. The solid lines represent two- 
dimensional runs, while the dashed lines correspond to their three- 
dimensional counterparts. The parameters of the two-dimensional 
simulations are described in Table Q] (runs 2D-7, 2D-8, and 2D-3 
for rrii/m e = 25, 100, and 400, respectively) and the ones of the 
three-dimensional runs are specified in Table [2] (runs 3D-1, 3D-2, 
and 3D-3 for rrii/m e = 25, 100, and 400, respectively). A con- 
vergence between the two- and three-dimensional results can be 
seen as mi/m e increases to more realistic values, implying a de- 
crease in the contribution of the "shock surfing" mechanism to the 
acceleration of electrons. 

and 2D-20 in Tabled]). Notice that in order to keep the 
same Ma = 7 the magnetic field also has to be reduced by 
a factor of 3.3. The resultant downstream electron spec- 
tra are shown in Figure [TOl for mi/m e = 25, 100, and 400 
(red, black, and green lines). These spectra essentially re- 
produce the ones corresponding to two-dimensional runs 
with v s h — 0.14c, shown in Figure |5] The same spectral 
index a dependence on m,/m e is obtained, with a de- 
crease in a when passing from mi/m e = 100 to 400. The 
overall normalization of the power law tails, however, 
seems to drop by a factor of ~ 3 with respect to the 
v s h = 0.14c cases. This result suggests that the mecha- 
nism for electron acceleration due to growth of whistler 
waves has a rather weak dependence on the shock veloc- 
ity, as long as Ma is kept constant. It is important to 
notice, however, that these measured spectral indices are 
for the same inclination angles that produce the lowest 
a in the case of v s h = 0.14c (Obu = 75°). Thus, in this 
case we are assuming that the angle Obu at which a is 
minimized does not depend on the shock velocity. This 
assumption needs to be confirmed by further exploration 
of the angle dependence of acceleration for different shock 
velocities. A thorough study aimed to clarify this point 
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Fig. 9. — The electric field along x, E x , in the shock transition re- 
gion of simulations 2D-7, 2D-8, and 2D-3, which only differ in their 
mi/rrie parameter given by mi/m e =25, 100, and 400, respectively. 
The field is shown at tw C)i = 10 (tw p , e = 2500, 5000, and 10000, 
respectively), and is normalized in terms of Eo, which corresponds 
to the convective electric field given by BoVi„/c (where Vi = 0.1c 
is the injection velocity of the upstream plasma as seen from the 
downstream frame). 

will be presented elsewhere. 

4.4. Alfvenic Mach Number Dependence 

We test the effect of varying Ma by running two- 
dimensional simulations that are analogous to the ones 
shown in Figure [SJ but also using the Mach numbers 
Ma — 3.5 and Ma = 14. We do this by changing 
the magnitude of the magnetic field while keeping the 
same shock velocity v s h = 0.14c. We use the mass ratios 
mi/m e — 25, 100, 400, and 1600, in order to simulta- 
neously take into account the mass ratio dependence for 
each of the studied Mach numbers. The parameters of 
the simulations are summarized in Table [TJ (where for 
m t /m e = 25, 100, 400, and 1600 the M A = 3.5 simu- 
lations are called respectively 2D-10, 2D-11, 2D-12, and 
2D-13, the M A = 7 simulations are called 2D-7, 2D-8, 
2D-3, and 2D-9, and the Ma = 14 simulations are called 
2D-14, 2D-15, 2D-16, and 2D-17). 

The corresponding downstream electron spectra for 
Ma = 3.5, 7, and 14 are shown in panel a), b), 
and c) of Figure [TTJ (with the red, black, green, and 
blue lines showing the spectra for rrii/m e — 25, 100, 
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Fig. 10. — The downstream electron spectra for the 2D-18, 2D-19, 
and 2D-20 simulations, with v sh = 0.042, Ma = 7, 8 B n = 75°, and 
mi/m e =25, 100, and 400 (depicted by red, black, and green lines, 
respectively). A dependence of the spectral index a on mi/m e 
similar to the one found for the v sh = 0.14 case is observed. 

400, and 1600, respectively). We see that, apart from 
the hardening of the spectra due to the increase in 
the mass ratio, the spectral index a also drops as the 
Mach number goes from 14 to 3.5. This tendency is 
consistent with our finding that the electron accelera- 
tion is driven by the growth of whistler waves in the 
foot of quasi-perpendicular shocks. Indeed, our analy- 
sis presented in t j3.ll showed that the growth of whistler 

waves requires MjJJ^mJmS^. < 1 dWu et al l 119831 ; 
IMatsukivo fc Scholer 1 120031 : ferasnoselskikh et al.ll2002D . 
In Figure [TTI we see that when this condition is not satis- 
fied (for instance, when Ma = 7 and m;/m e = 25, shown 
by the red line in plot b) we obtain a rather soft spec- 
trum, with a ~ 6. The acceleration efficiencies, in terms 
of the fraction of non-thermal particles, follow a similar 
trend, with the number of non-thermal particles being 
larger for smaller values of Ma/ '(mi/me) 1 / 2 . The typical 
fraction of non-thermal particles in our runs ranges be- 
tween ~ 2 to 10%, with corresponding energy fractions 
of ~ 10 to 30%EJ. 

For each Alfvenic Mach number, we explored the angle 
dependence of the electron acceleration by using several 
values for Obu- Thus, the spectra presented in Figure [Til 
correspond to the hardest spectra for each Ma, whose an- 
gles are: 9 Bn = 75° for M A = 3.5, 6 Bn = 75° for M A = 7, 
and 0Bn — 60 for Ma — 14. However, these angles were 
determined using a single mass ratio irii/m e = 400. Thus 
our results assume that the angle of the hardest spectral 
index depends only weakly on the used mass ratio. This 
assumption needs to be confirmed by further exploration 
of the angle dependence of this acceleration mechanism 
for different values of mj/m e . Despite this uncertainty, 
our results show, in general, a hardening of the electron 
spectra as the Ma/ (mi/me) 1 / 2 ratio decreases. This can 
also be seen from Figure[HJ which shows a compilation of 
the spectral indices for the cases depicted in Figure [TTI as 
a function of the mass ratio and Mach number. The yel- 
low dashed line, which corresponds to Ma = (mi/me) 1 / 2 , 
clearly separates the regions of small a from regions of 

4 These percentages are estimated by considering the electrons 
with energies in the range where the power-law component of the 
spectrum is larger than the thermal component. 
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Fig. 11. — The downstream electron spectra are shown for simulations with M A = 3.5, 7, and 14, in panels a), 6), and c), respectively. 
The variations in M A are obtained by only changing the magnitude of the magnetic field, so in all the cases the shock velocity has the 
value v s h = 0.14c. The realistic mass ratio mi/m e = 1600 was also included, so the red, black, green, and blue lines represent cases with 
mi/m e = 25, 100, 400, and 1600, respectively. We can see from these spectra the overall trend to have smaller spectral index a either when 
the Mach number is reduced or when the mass ratio is increased. The parameters of each of these simulations are compiled in Table [T] 
(where for mi/m e = 25, 100, 400, and 1600 the M A = 3.5 simulations are called respectively 2D-10, 2D-11, 2D-12, and 2D-13, the M A = 7 
simulations are called 2D-7, 2D-8, 2D-3, and 2D-9, and the M A = 14 simulations are called 2D-14, 2D-15, 2D-16, and 2D-17). 
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Fig. 12. — A compilatio n o f the spectral indices for the 12 simu- 
lations depicted in Figure [TT] shown as a function M A and mi/m e . 
The yellow dots correspond to the actual grid point locations, 
while the rest of the diagram is colored using linear interpola- 
tion of the measured values of a. The yellow dashed line cor- 
responds to M A = (mi/me) 1 / 2 . We see that regions of hard 
and soft spectra tend to be separated by this line, which marks 
the theoretical estimat e for the limit between growth and suppres- 
sion of whistler waves HWu et aLlflMllMatsukivo fc Scholer |[2003l ; 
IKrasnoselskikh et al.1120021) . 

rather soft spectra. In the cases with realistic mass ratio 
(m,/m e = 1600) the spectral index goes from a = 2.7 to 
4 as Ma changes from 3.5 to 14. We will discuss the con- 
sequences that this result has for electron acceleration in 
SNR shocks below. 

Finally, we emphasize that this acceleration mecha- 
nism does not show a significant dependence on the cho- 
sen value of (3 e . This is verified by directly comparing run 
2D-3 and 2D-21, which only differ in their j3 e (/3 e =0.5 and 
0.05, respectively) and show essentially no difference in 
their final electron spectra. 

5. SUMMARY AND CONCLUSIONS 

In this paper we studied the problem of electron accel- 
eration in non-relativistic electron-ion shocks, using two- 



and three-dimensional PIC simulations. By systemati- 
cally exploring the space of shock parameters, and using 
realistic ion to electron mass ratios, we identify a new 
acceleration mechanism, that is able to accelerate elec- 
trons starting from fairly low temperatures. Thus, it 
constitutes a possible candidate to solve the well known 
electron "injection problem" of the diffusive shock accel- 
eration (DSA) theory. This mechanism confirms the idea 
that non-thermal electron acceleration can be caused by 
electron scattering due to foot waves, which are produced 
by the electron-ion counter-str eaming (see, for example, 
ICargill fc Papadopoulos 1 0.9881 ) . We have described the 
physics of this mechanism, and identified the physical 
regime where this acceleration is most efficient. The 
acceleration occurs preferentially in quasi-perpendicular 
shocks, and is driven by oblique whistler waves excited 
in the shock foot. 

We found that, for simulations that only differ in their 
ion to electron mass ratio mj/m e , the spectral index 
tends to harden as the m,-/m e grows from 25 to 1600. 
This trend is confirmed for Ma — 3.5, 7, and 14, and 
for shock velocities v s h — 0.14c and 0.042c. On the 
other hand, simulations that only differ in their Alfvenic 
Mach numbers also show a progressive hardening of their 
spectra as Ma is reduced from 14 to 3.5. This mass 
ratio/ Alfvenic Mach number dependence of the acceler- 
ation is consistent with theoretical arguments suggest- 
ing that the growth of whistler waves in the foot of 
quasi-perpendicular shocks would be favored when Ma < 
(m.i/me) 1 / 2 . Although the physics of this whistler exci- 
tation is still a subject of debate, this condition appears 
to hold independe ntly of whether the whistlers are gener- 
ated by the MTSI (|Wu et al.lll983HMatsukivo fc Scholer I 
120031 l2006h . or if they are just explained as an intrinsic 
component of the structure of quasi-perpendicular shocks 
(Krasnoselskikh et al. 2002). We also found a strong de- 
pendence of the acceleration on the angle between the 
magnetic field and the shock normal, 0Bm which needs 
to satisfy 8b u ^ 90°. O n the other hand, we found that 
the shape of the spectra do not depend significantly on 
the shock velocity, although a factor of ~ 3 decrease in 
the normalization of the non-thermal part of the distri- 
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TABLE 1 

Parameters for the two-dimensional simulations 



Run 


c/0J P! e 




(3e = Pi 




M A 


e Bn 


m,i/m e 




2D-1 


10 


102 


0.5 


0.14 


7 


90° 


400 


100 


2D-2 


10 


102 


0.5 


0.14 


7 


80° 


100 


100 


2D-3 


10 


102 


0.5 


0.14 


7 


75° 


400 


100 


2D-4 


10 


102 


0.5 


0.14 


7 


70° 


400 


100 


2D-5 


10 


102 


0.5 


0.14 


7 


60° 


400 


100 


2D-6 


10 


102 


0.5 


0.14 


7 


45° 


400 


100 


2D-7 


15 


34 


0.5 


0.14 


7 


75° 


25 


100 


2D-8 


10 


77 


0.5 


0.14 


7 


75° 


100 


100 


2D-9 


10 


102 


0.05 


0.14 


7 


75° 


1600 


10 


2D-10 


1 


102 


n ni 


14 


o.o 


( u 




1 


2D- 11 


10 


102 


0.01 


0.14 


3.5 


75° 


100 


10 


2D-12 


10 


102 


0.01 


0.14 


3.5 


75° 


400 


10 


2D-13 


10 


102 


0.01 


0.14 


3.5 


75° 


1600 


10 


2D-14 


10 


102 


0.2 


0.14 


Id 


60° 


25 


10 


2D-15 


10 


102 


0.2 


0.14 


14 


60° 


100 


10 


2D-16 


10 


102 


0.2 


0.14 


14 


60° 


400 


10 


2D-17 


10 


102 


0.2 


0.14 


14 


60° 


1600 


10 


2D-18 


15 


34 


0.005 


0.042 


7 


75° 


25 


100 


2D-19 


15 


51 


0.005 


0.042 


7 


75° 


100 


32 


2D-20 


10 


77 


0.005 


0.042 


7 


75° 


400 


65 


2D-21 


10 


102 


0.05 


0.14 


7 


75° 


400 


100 


Note. 


— We list the electron 


skin depth 


c/w PiC in 


terms of number of grid 


cells, 



the transverse size of the simulation box in terms of c/cj P;e , the beta parameter of 
the different plasma particles (3j (= pj / Bq /87r, where pj is the pressure of particle 
"j"), the upstream medium shock velocity, v s h, the Alfvenic Mach number, Ma, the 
angle between the upstream magnetic field and the shock normal, Qun, the ion to 
electron mass ratio, mi/m e , and the total number of particles per cell, N ppci in the 
simulation. 



TABLE 2 

Parameters for the three-dimensional simulations 



Run c/ui Pi e L y X L z /(c/cjp, e ) 2 f) e = Pi v sh /c Ma 9 Bn rrii/rn e N ppc 

3D-1 5 26 x 3 0.5 0.14 7 75° 25 100 

3D-2 5 77 x 3 0.5 0.14 7 75° 100 25 

3D-3 5 51 x 3 0.5 0.14 7 75° 400 50 

Note. — We list the electron skin depth c/uj p ^ e in terms of number of grid cells, the 
transverse size of the simulation box in terms of c/u PiE , the beta parameter of the different 
plasma particles fij (= Pj/Bq/St:, where pj is the pressure of particle "j"), the upstream 
medium shock velocity, v 8 h, the Alfvenic Mach number, M^, the angle between the upstream 
magnetic field and the shock normal, 0Bm the ion to electron mass ratio, m,/m e , and the 
total number of particles per cell, N ppc , in the simulation. 



bution was observed when the shock velocity was reduced 
from v 8 h = 0.14c to v s h = 0.042c. 

This dependence of the acceleration on the shock pa- 
rameters can be explained in terms of two requirements: 
1) the shock needs to be able to excite whistler waves 
in the foot, which translates into the condition Ma < 
(mj/me) 1 / 2 , and 2) the accelerated electrons need to 
stay in the shock foot for a time long enough to increase 
their chances of being scattered by the whistler waves. 
The second requirement explains why the efficiency of 
the acceleration drops substantially when Obu = 90°. 
In that case, energized electrons cannot move along the 
shock normal, n, since they are tied to the magnetic 
field. Thus, they are rapidly advected into the down- 
stream medium of the shock. This does not happen in 
an oblique shock, where accelerated particles can develop 
a significant velocity along h. Given that the growth 
of whistler waves does not depend on the shock veloc- 
ity, v s h, the weak dependence of the acceleration on v s h 
implies that the ability of the accelerated electrons to 
stay in the shock foot is also independent of v s h- We 



can verify this point as follows. The initial electron en- 
ergization is driven by the electric field of the whistler 
waves (in particular, its component parallel to the mag- 
netic field). This electric field has a typical length scale 
of ~ 10c/w Pi e and a maximum amplitude of ~ lO-Eo 
(where Eq = Bov s h/c). Thus, the energy gain in each 
scattering will be given by ~ 100 eE c/uj Pte . There- 
fore, an electron that moves along Bq with that en- 
ergy, will have x— velocity larger than the shock veloc- 
ity if 100 eEo c/u!p te cos(dBn) 2 m eV 2 h' which implies 
100cos(6>b„) 2 > MA/imi/rrie) 1 / 2 . We can see that this 
condition is independent of the shock velocity, and con- 
firms that Qbu and the M J 4/(mi/m e ) 1 / 2 ratio are the 
crucial parameters for the acceleration mechanism. 

The maximum electron energy that we measure is con- 
sistent with the electron Larmor radii Rl, £ being com- 
parable to that of the ions i?L,i, as can be verified from 
Figure 1 1 31 where we plot the downstream spectrum of 
electrons for run 2D-3 (M A = 7, v sh = 0.14c, 9 Bn = 75°, 
and mi/m e = 400). The electron spectrum is shown as 
black line, and the thermal and power law tail fits are 
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Fig. 13. — The downstream energy spectra for electrons and ions 
(black and green lines, respectively) are shown for the shock tran- 
sition region of run 2D-3 (Ma = 7, v a h = 0.14c, dsn = 75°, and 
mi/m e = 400). The thermal and power law tail fits are shown 
in red dashed lines. The vertical dotted line marks the electron 
Lorentz factor corresponding to equal electron and ion Larmor radii 
(r e = (mi/m e ){v sh /c)). 

depicted in red dashed lines (a « 3.6). The vertical dot- 
ted line marks the electron Lorentz factor corresponding 
to R.L,e = Rl,zi which coincides with the maximum en- 
ergy of the electrons. However, this maximum energy 
is statistically limited by the number of macroparticlcs 
used in the simulations. Thus, higher energy electrons 
may in principle be possible. In Figure [13] we also de- 
pict the corresponding downstream ion spectrum, which 
shows incomplete thermalization and the absence of a 
power law tail. 

We also explored the contribution to acceleration given 
by the previously proposed "shock surfing" mechanism, 
which is driven by the excit ation of Buneman waves i n 
the leading edge of the foot ([Amano fc Hoshino 1 12009). 
We find that its contribution is strongly dependent on 
the artificial mass ratio mi/m e , becoming negligible for 
realistic values of rrii / m e . 

From the observational viewpoint, our results are con- 
sistent with in-situ electron spectrum measurements in 
the Earth's bow shock, in which t ypically Ma ~ 5 — 10. 
For instance. iGosling et al.l (|1989T ) reported on spectrum 
measurements performed on the International Sun-Earth 
Explorer (ISEE) 1 and 2 spacecraft, showing that non- 
thermal electrons are commonly found in the down- 
stream medium of the quasi-perpendicular portions of 
the shock, but are rarely observed in the quasi-parallel 
portions. They also found typical spectral indices of 
a ss 3 — 4, with the power-law part of the spectrum 
extending smoothly out of the thermal part of the dis- 
tribution, which is essentially what we observe in our 
si mulations. Similar results were more recently reported 
bv lOka et all (2006), using data from the Geotail space- 
craft. We emphasize, however, that the application of our 
results to the Earth's bow shock measurements is based 
on simulations with v s h = 0.14c and 0.042c, which are a 
factor of at least ^100 larger than the shock velocities 
typically found in the solar system. This extrapolation is 
made due to the finding that the crucial parameters for 
the acceleration due to whistlers are Ma and Qbu, with 
the shock velocity dependence being much weaker. This 
extrapolation remains to be confirmed by shock simula- 



tions with orders of magnitude smaller v s h (or va)- 

In the context of electron acceleration in SNR shocks, 
we can confirm the potential importance of this injec- 
tion mechanism by estimating the Alfvenic Mach num- 
ber necessary to explain the typical fraction of parti- 
cles injected into the DSA process, T)i n j. According to 
the modeling of broadband observations of SNRs, rj in j 
ranges between ~ 10 -4 , for magnetic fields amplified 
only by shock compression (B ~ 17fj,G), and ~ 10~ 6 , 
for significant magnetic amplification consistent with re- 
cent X-ray observations (B ~ 13 0/j.G 1 ) (see, for example, 
Zirakashvili & Aharonian 2010). Also, in order to be 
injected, the particles need to have Larmor radii close 
to the size of the shock transition region, whose char- 
acteristic scale corresponds to the Larmor radii of the 
ions. Thus, we can estimate the maximum slope of the 
non-thermal part of the spectrum necessary to satisfy 
these requirements. If we assume that the normalization 
of the non-thermal spectrum is such that the power law 
tail dominates for energies very close to the peak of the 
thermal distribution (which is approximately what we 
see in our simulations), we get that the necessary spec- 
tral index should satisfy 

awl- log(r]i nj )/log(c/v sh ). (1) 

Thus, for a typical shock velocity of v s h = 3000 km/s, 
and f]i n j — 10~ 6 , we get that a « 4. From the results 
summarized in Figure [T2l we sec that, for realistic mass 
ratios (mj/m e = 1600), the maximum Mach number that 
would give a < 4 corresponds to Ma = 14. This estimate 
shows that the injection mechanism presented in this pa- 
per is a viable solution to the injection of electrons into 
the DSA mechanism in SNR shocks only if the Alfvenic 
Mach number is smaller than ~ 20. This result implies 
that significant magnetic field amplification must occur 
in the upstream region of SNR shocks. Indeed, if the up- 
stream field is not amplified, the typical Alfvenic Mach 
number for SNR shocks would be ~ 300, for a shock with 
v s h = 3000 km/s propagating in a rii = 1 cm~ 3 plasma 
with a 3/xG ISM field. This is interesting considering that 
significant magnetic amplification in SNR shocks based 
on X - ray observations has be en recently reported {Ballet 
l2006t lUchivama et alj [2007) . These observations show 
the existence of thin, non-thermal rims, which are inter- 
preted as synchrotron emission by TeV electrons acceler- 
ated at the shocks. The rapid variability and thinness of 
the rims (which depend on the synchrotron cooling time 
of the electrons) have allowed to estimate the strength 
of the field, suggesting downstream amplitudes ~ 100 
times larger than typically expected in the ISM of the 
Galaxy. This implies an amplification factor of the up- 
stream field of ~ 25, considering that the field compres- 
sion at the shock itself would contribute an extra factor 
of ~ 4 to the total downstream field growth, making it 
possible to have Ma < 20 in these shocks. Also, there 
are theoretical reasons to believe that this field growth 
could happen in the upstream medium of the shocks, due 
to streaming instabilities driven by cosmic ra ys (CRs) as 
they get accelerated in these environmen ts (jBelll I2004L 
2005: IRiauelme fc Spitkovskv1[200l l2010h . 

The Ma % 20 condition for electron injection is also 
supported by a recent PIC study of a two-dimensional, 
perpendicular shock with in-plane magnetic field and pa- 
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rameters: Ma ~ 130, Vi n = 0.25c, and mi/m e = 30 
(|Kato fe Takabe Il2010l ). This study showed the absence 
of electron acceleration and whistler waves in the foot of 
the shock. Indeed, the shock foot is dominated by the 
Weibel instability, which plays a fundamental role in the 
formation of the shock itself. Although the strict per- 
pendicularity of the shock (9bu = 90°) may contribute 
to reducing the electron acceleration, we believe that the 
suppression of the whistler waves in this low magnetiza- 
tion shock would make electron injection unlikely for any 
value of Obu- 

Thus, this result reinforces the idea of very large mag- 
netic amplification associated with the synchrotron rims 
seen by the X-ray observations of SNRs. It is important 
to emphasize that in this paper we have concentrated 
only on the injection part of the acceleration process. Af- 
ter the injected electrons reach energies such that they 
can move diffusively in the shock vicinity, they would be 
further accelerated via the DSA up until the ~ 10 — 100 
TeV energies inferred from observations. The DSA part 
of the acceleration process is not captured in our simu- 
lations due to the lack of upstream turbulence. 

As seen from Figure 1131 ions are not accelerated in the 
quasi-perpendicular configurations studied in this work. 
On the other hand, in quasi-parallel shocks, Alfven waves 
would interact resonantly with returning ions, providing 
the pitch angle scattering necessary to confine them to 
the shock vicinity, and all owing their acceleration v ia the 
first order Fermi process (|Kulsrud fe Pearcd U969L Gar- 
gate & Spitkovsky 2010, in prep.). Thus, given that the 
acceleration of electrons and ions would require differ- 
ent magnetic orientations, in principle, the acceleration 
of these two species would not be possible in the same 
regions of the shocks. This fact poses a consistency prob- 
lem when applying the electron injection mechanism pre- 
sented here to the case of SNR shocks, given that in SNR 
conditions the acceleration of electron requires the pres- 
ence of relativistic ions to amplify the held. This problem 
would be solved if the magnetic growth could also involve 
magnetic field reorientation. This is actually what has 
been found by numerical studies of the CR streaming 
instabilities, which are characterized by very turb ulent 
magnetic configurations in their non-linear sta te (|Belll 
12001 [20051 fRiauelme fc Spitkovsky 1 120091 l2Q10h . Also, 
these instabilities would produce magnetic field fluctu- 
ations on scales comparable to the Larmor radii of the 
CRs, R,L,cr- Therefore, given that the electron injec- 
tion mechanism presented here operates on length scales 
comparable to a few c/uj Pt i (which is much smaller than 
RL,cr)i h is plausible that electron acceleration can hap- 
pen locally, in regions where B is quasi-perpendicular to 
the shock normal. 

The likely global picture of the shock acceleration pro- 
cess then unfolds as follows. On large scales upstream 
of the efficiently accelerating shock, the magnetic field 
is likely quasi-parallel to the shock normal. This con- 
figuration is generally conducive to ion (cosmic ray) ac- 
celeration and subsequent escape. Escaping cosmic rays 
amplify and reorient the upstream magnetic field closer 
to the shock via current-driven instabilities. The re- 
sulting magnetic turbulence advected to the shock will 
be roughly isotropic, with regions of quasi-perpendicular 
magnetic field intermittently crossing the shock. This lo- 



cally transverse field will be efficient in injection of elec- 
trons via the whistler mechanism. Once pre-accelerated, 
the electrons will join the ions in DSA on scales larger 
than the shock foot in the turbulence that is driven by 
the cosmic rays. This speculative picture underscores 
the interdependence of electron and ion acceleration, as 
without cosmic rays the amplification and reorientation 
of magnetic field needed for electron injection would be 
hard to achieve. Observationally, this scenario is con- 
sistent with the large scale magnetic field in SN1006 
pointing along the axis that connects the "polar caps" 
of bright non-thermal emission (Ballet 2006). Also, the 
degree of linear polarization should be smaller in SNR 
rims with strong synchrotron emission, as we expect the 
field direction near t he shock to be random ized by the 
amplified turbulence (jStroman fc Poh l 2009). On larger 
scales, the radial magnetic fields inferred from radio po- 
larization measurements of synchrotron-emitting regions 
of SNR shocks are also consistent with our model (Dickel 
et al. 1991, DeLaney et al. 2002). 

The results presented in this paper correspond to only 
a partial exploration of the space of shock parameters. 
Studying this problem using realistic mass ratios is com- 
putationally expensive, so testing all the possible com- 
binations of parameters of interest is complicated and 
will require further work. For instance, one of our basic 
assumptions is that the angle Osn at which the acceler- 
ation is maximized has a weak dependence on the ion 
to electron mass ratio. We believe that this assumption 
needs to be investigated by performing a more complete 
exploration of the acceleration efficiency dependence on 
0Bn and rrii/m e . In a similar way, the dependence of 
the acceleration efficiency on the shock velocity v s h re- 
quires further study. We found that, when passing from 
v s h = 0.14 to 0.42, the slope of the power law part of the 
electron spectra for the mass ratios mi/m e = 25, 100, 
and 400 does not change significantly. However, this re- 
sult was obtained only in the case of Osn — 75° and 
Ma = 7. We think that it would be interesting to con- 
firm this weak velocity dependence by performing a thor- 
ough exploration of the shock parameter space at lower 
velocities. Finally, comparing two-dimensional simula- 
tions with analogous three-dimensional runs where the 
typical whistler length scale (~ 10c/uj Pje ) is resolved in 
the third dimension would also be useful to confirm our 
results and to enquire about possible additional three- 
dimensional effects. Although we do not expect the qual- 
itative picture presented here to change significantly, we 
believe that these studies would give us more accuracy 
in our estimates of the electron injection efficiencies as a 
function of the shock parameters. 

With these caveats in mind, in this paper we have 
identified a new possible mechanism for electron injec- 
tion into the DSA process in non-relativistic shocks. 
The obtained spectra are consistent with in-situ mea- 
surements of electron energy distributions at the Earth's 
bow shock. Also, our results would explain the observed 
electron acceleration in SNR shocks, implying very large 
upstream magnetic field amplification in these environ- 
ments. Thus, if this mechanism proves to be the only 
possible solution for electron injection in non-relativistic 
blast waves, it would reinforce the inferred strong con- 
nection between particle acceleration and magnetic field 
amplification in SNR shocks. 
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Fig. 14. — The downstream electron spectra at tui c i = 10 for the two-dimensional simulations 2Db-l (red line), 2Db-2 (black line), and 
2Db-3 (green line), described in Table|3] The simulations are characterized by having a magnetic field quasi-perpendicular to the simulation 
plane, and only differ in their mass ratios, which have values mi/m e = 25, 100, and 400, respectively. 
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APPENDIX 

MASS RATIO DEPENDENCE OF THE "SHOCK SURFING" ACCELERATION 

In this appendix we exp lore the mass ratio dependence of the electron energization due to shock surfing acceleration 
ijAmano &: Hos hino 2009) . This acceleration is driven by the presence of electrostatic waves, produced by the Buneman 
instability (|Buneman| [l958). in the foot of quasi-perpendicular shocks. The Buneman waves grow in the leading edge 
of the foot due to the relative velocity between the upstream electrons and the shock-reflected ions, and have a typical 
scale comparable to ~ c/w Pje . The small length scale of these waves, which can be comparable or even smaller 
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Fig. 15. — The electric field along x, E x , in the shock transition region of simulations 2Db-l (panel a), 2Db-2 (panel b), and 2Db-3 (panel 
c), which are described in Table [3] These simulations are characterized by having a magnetic field quasi-perpendicular to the simulation 
plane, and only differ in their mi/m e parameter given by m;/m e =25, 100, and 400, respectively. The field is shown at tui c , = 10, and is 
normalized in terms of Eq, which corresponds to the convective electric field given by BoVi n /c (where = 0.1c is the injection velocity of 
the upstream plasma as seen from the downstream frame). 



TABLE 3 

Parameters for the two-dimensional simulations with magnetic field 
quasi-perpendicular to simulation plane 



Run 




Ly/c/u)p t< . 


Pe = A 


Vsh/c 


M A 




rm/m e 




2D-lb 


15 


26 


0.5 


0.14 


7 


75° 


25 


100 


2D-2b 


15 


34 


0.5 


0.14 


7 


75° 


100 


100 


2D-3b 


10 


51 


0.5 


0.14 


7 


75° 


400 


100 


2D-4b 


10 


51 


1.5 


0.14 


21 


75° 


25 


100 


2D-5b 


10 


26 


4.5 


0.14 


21 


75° 


100 


100 


2D-6b 


15 


34 


0.045 


0.042 


7 


75° 


25 


100 


2D- 7b 


10 


26 


0.045 


0.042 


7 


75° 


100 


100 


Note. 


— We list the electron 


skin depth 


c/w Pifi in 


terms of number of grid 


cells, 



the transverse size of the simulation box in terms of c/u P;e , the beta parameter of 
the different plasma particles 0j (= pj / Bq /8-7T , where pj is the pressure of particle 
"j"), the upstream medium shock velocity, v s h, the Alfvcnic Mach number, Ma, 
the angle between the upstream magnetic field and the shock normal, Osm the ion 
to electron mass ratio, m,/m e , and the total number of particles per cell, N ppc , in 
the simulation. 

than the electron Larmor radii, makes the Buneman waves a potentially important means for electron scattering and 
energization in the foot of quasi-perpendicular shocks. As the upstream electrons encounter these waves, they gain 
energy due to the work performed by the fluctuating electric field, combined with the convective field of the upstream 
plasma. The wave vector of the fastest growing mode is parallel to the velocity of the beam of returning ions, so these 
waves are best studied either by three-dimensional simulations or by two-dimensional runs where the magnetic field is 
quasi-perpendicular to the simulation plane, so that the gyrational motion of the ions is resolved on the plan^B 

Using a two-dimensional simulation with shock param eters Vj n = 0.2c (v s h ~ 0.3 c), Ma = 14, 9bu = 90° (with Bq 
perpendicular to the simulation plane), and mi/m e = 25. IAmano fc Hoshino I (|2009f ) showed that an electron spectrum 
with index a ~ 2-2.5 can be produced in the shock transition region due to this mechanism. 

In this appendix we show, however, that the importance of this shock surfing acceleration decreases as the mass ratio 
increases to more realistic values. As an example, let us compare the downstream spectra of three simulations with 
the same shock parameters Vi n = 0.1c (v s h w 0.14c), Ma = 7, Osn = 75° (and the angle between Bq and z equal to 
15°), but different mass ratios: mi/m e = 25, 100, and 400 (the rest of the parameters are specified in Tabled]). Figure 
1141 shows the spectra at tcj c ,i = 10 measured in the downstream medium. Although in all three cases a non-thermal 
component can be seen, the maximum energy reached by the electrons is about the same (r e - 1 a 1), despite the 
factor of 16 difference in the used mass ratios. In particular, in the rrii/m e — 400 case, the closeness between the 
peak of the thermal part of the distribution and the maximum energy of the accelerated particles reduces significantly 
the relative fraction of non-thermal particles compared with rn,i/m e = 25 and 100 cases. Also, notice that in these 

5 As in the main part of the paper, in this appendix the two- 
dimensional simulations will be in the xy plane. 
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two-dimensional simulations the acceleration due to whistler waves gets suppressed. This is because these simulations 
only allow the wave vector of the whistlers to be quasi-perpendicular to the magnetic field, inhibiting their growth and 
their effect on the electron energization. 

In terms of the evolution of the Buneman waves, we also observe a decrease in their relative amplitude. This can 
be seen from Figure [T5l where the electric field along the a;— axis is depicted for the same simulations shown in Figure 
1141 The Buneman waves correspond to the fluctuations on ~ c/w Pie scale that appear in the leading edge of the shock 
(at x/c/ui Pye = 115 — 120, 240-245, and 365-370, in panels a), b), and c), respectively). The black dashed line in the 
three panels marks the position of the shock density peak (overshoot) in the three simulations. When compared to 
the convective electric field (E = B v in /c), we see that the maximum electric field of the Buneman waves is reduced 
by a factor of ~ 4 when mi/m e passes from 25 to 400. 

The same maximum energy for the accelerated electrons is observed in similar simulations but using Ma = 21 
(simulations 2Db-4 and 2Db-5). Also, when trying cases similar to the ones shown in Figures IT4l and [TBI but using 
v s h — 0.042 (simulations 2Db-6 and 2Db-7), a decrease in the maximum energy to T e — 1 » 0.1 is seen, with no change 
by varying rrii/m e . 

Thus, in general, the maximum energy attainable by the shock surfing mechanism appears to depend mainly on the 
shock velocity, and shows no changes for different values of m,/m e . Therefore, as the mass ratio mi/m e is increased, 
the peak of the thermal part of the distribution gets closer to this maximum energy. This implies that the fraction 
of non-thermal electrons is reduced when rrii/m e approaches more realistic values. We conclude that this mechanism 
would not contribute significantly to electron acceleration in the fully realistic case of mi/m e = 1836. 



